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Estimating column density from ammonia (1,1) emission in 
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ABSTRACT 

We present a new, approximate method of calculating the column density of ammonia in 
mapping observations of the 23 GHz inversion lines. The temperature regime typically found 
in star forming regions allows for the assumption of a slowly varying partition function for 
ammonia. It is therefore possible to determine the column density using only the (J=1,K=1) 
inversion transition rather than the typical combination of the (1,1) and (2,2) transitions, with 
additional uncertainties comparable to or less than typical observational error The proposed 
method allows column density and mass estimates to be extended into areas of lower signal to 
noise ratio. We show examples of column density maps around a number of cores in the W3 
and Perseus star-forming regions made using this approximation, along with a comparison to 
the corresponding results obtained using the full two-transition approach. We suggest that this 
method is a useful tool in studying the distribution of mass around YSOs, particularly in the 
outskirts of the protostellar envelope where the (2,2) ammonia Une is often undetectable on 
the short timescales necessary for large area mapping. 
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1 INTRODUCTION 

The usefulness of ammonia as a tracer of conditions in the interstel- 
lar medium (ISM) has been recognised for many years (e.g. ???). 
Ammonia is one of the most abundant molecular species in inter- 
stellar clouds, with typical observed molecular abundance ratios 
ranging from 10"' - 10"' (e.g. ???). The excitation of the (1,1) ro- 
tational inversion transition of ammonia at 1.3 cm at typical molec- 
ular cloud temperatures (10 - 100 K) is overwhelmingly due to col- 
lisional processes (?) and thus, above a critical density of ~4x 10^ 
cm"^ (?), is sensitive to the temperature of the region under ob- 
servation. These densities are typical for star-forming regions (e.g. 
?). 

The hyperfine structure of the NHj (1,1) transition allows the 
optical depth associated with the line to be determined from a sin- 
gle observation, while other molecular species require separate ob- 
servations of an (often rare) isotopologue. Observations of the (2,2) 
transition also allow measurement of the rotational and kinetic tem- 
perature of the gas through the line ratios. The proximity of the 
(1,1) and (2,2) transitions in frequency means that the two can of- 
ten be measured in the same instrumental bandpass (i.e. in a single 
observation), eliminating or greatly decreasing uncertainties arising 
from calibration and pointing eiTors. 

The combination of optical depth and temperature measure- 
ments allow the column density of molecular gas to be determined 
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through commonly used methods (?) and the assumption of LTE. 
Through this application, the distribution of molecular gas column 
density in the dense ISM may be mapped. Such mapping simulta- 
neously determines the velocity dispersion of the region under in- 
vestigation, providing virtually all of the physical information that 
may be determined in a single set of observations. For these reasons 
ammonia has become one of the most useful observational tools in 
modem studies of young stellar objects (YSOs), (e.g., ???). The 
ramifications for the study of star formation include the determina- 
tion of virial ratios (c.f. ?), the relative abundances of nitrogen- vs. 
carbon-bearing molecular species (e.g. ?) and the relationship be- 
tween the distribution of submillimetre continuum emission from 
dust to the emission from the associated molecular gas (???)). 

Recent developments in instrumentation, such as the K-band 
Focal Plane Array (KFPA) at the Green Bank Telescope (GBT) in 
West Virginia, mean that ammonia emission may now be mapped 
over larger areas than were previously possible within reasonable 
time constraints. However, the maps resulting from such observa- 
tions are likely to suffer from poor signal-to-noise ratios (SNRs) 
in comparison to pointed observations, away from the position of 
peak emission. The reasons for this are that ammonia observations 
of star-forming regions are often based upon previous submillime- 
tre continuum emission maps and pointed observations are usu- 
ally made towards suspected high column-density lines of sight, 
while mapping observations are probing fainter regions. Also, map- 
ping observations are often time-consuming in comparison to sin- 
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Figure 1. Example NH3 (1,1) (bottom) and (2,2) (top) spectra taken towards 
W3 at the coordinates 02:25:33.0 +61:13:21 (J2000) 



gle pointings and so shorter integration times per-pixel are often 
required. 

Recent results from large scale surveys such as the Red MSX 
Source (RMS) Survey and the Herschel infrared Galactic (HiGAL) 
Plane Survey indicate consistent average temperatures over a wide 
range of star-forming environments of ~ 16 - 25 K (????). Over 
this temperature range the NH3 (2,2) transition (corrected) antenna 
temperature {Tl^^^i)) comparable to or greater than the line tem- 
peratures of the (1,1) hyperfine quadrupole transitions (often called 
'satellite' lines, labelled 'A', 'B', 'D' and 'E' in Fig.|T]l. The outer 
('A' and 'E') and inner ('B' and 'D') satellite quadrupoles have in- 
tensity ratios to the main quadrupole ('C') of 0.22 and 0.28 respec- 
tively (?). This corresponds to ratios of peak antenna temperature of 

~0.3 and, for a main line optical depth (Tmd n) of 1, ranges 

^Am(l,]) 

from 0.3 to 0.8. However, below the same temperature range of 16 

r* 

- 25 K, falls below 0.3. Therefore, determination of optical 

' Am(l,l) 

depth, which requires detection of the (1,1) hyperfine quadrupole 
transitions, is likely to be ascertained over a greater spatial region 
than the rotational/kinetic temperature, which requires detection of 
the (2,2) main quadrupole. This limits the area over which column 
density may be determined to that in which the weaker (2,2) line is 
detected, despite more extensive detection of ammonia in the (1,1) 
transition. 

In this paper we propose an approximate method of estimating 
column density along lines of sight in which only the (1,1) transi- 
tion is detected. The validity of the method requires only that the 
uncertainty introduced by variation of derived column density with 
kinetic temperature is less than typical observational errors over the 
range of temperatures found in star-forming regions. 



2 METHOD 

We calculate the column density of ammonia molecules in the (1,1) 
transition following equation A4 of ?, 



A'(l,l)^ 
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where vq, Tex and t are the line rest frequency, measured 



excitation temperature and optical depth respectively. A is the 
Einstein spontaneous emission coefficient and h and ks are the 
Planck and Boltzmann constants respectively. Here jT(v)dv = 

T,oi, where cr v is the measured line (velocity) disper- 
sion and T,ot is the total measured line optical depth, both deter- 
mined through line fitting. 

The total column density is determined by scaling the column 
density in the (1,1) state by the ratio of the overall partition func- 
tion to the individual (1,1) partition function (Z/Z(J=1)), i.e. the 
probability of occupation of the (1,1) state. 

The partition function for the metastable (J = K) states of the 
ammonia molecule is 



Z = (2J+ l)S(J)exp- 



-h [BJ(J+ l) + iC-B)V] 

knTk 



(2) 



Here S(J) is taken as 2 for J=3, 6, 9 ... and 1 otherwise, ac- 
counting for the extra statistical weight of ortho-NHs over para- 
NH3. The other terms in the partition function are, B and C - 
the rotational constants of the ammonia molecule, 298117.06 and 
186726.36 MHz, respectively (?); T^^ is the kinetic temperature. 

It can be seen then that the total column density of ammo- 
nia is dependent upon kinetic temperature (T^) in two ways, firstly 
through the contribution to the partition function and then in a more 
circuitous fashion as follows. In the assumption of LTE we have 
= T];, therefore, the population of the (1,1) state is also depen- 
dent upon T];. This is incorporated into our calculation of column 
density in equation [T] in using as the temperature defining the 
relative population of states (from the Boltzmann distribution func- 
tion). It should be noted that measurements of Tex are often lower 
than T^^, probably due to clumping of the gas on scales smaller than 
the telescope beam and the fact that T^.^ is derived from the single 
(1,1) transition. This inequality means that it is usually not valid 
to use Tex as an estimate of r^, particularly in single-dish observa- 
tions. 



2.1 Temperature Range 

Kinetic temperature is related to the rotation temperature through 
the following equation established by ? and updated with revised 
coUisional rates by ? 



l + iln[l+0.6exp(-J|^)]' 



(3) 



Where To = ~ 41. 5K. It should be noted that, for Tt, > 

Tq, r,. may be significantly underestimated (?). 



Tf is itself determined through the equation 
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(?), where T^^ represents the measured (corrected) antenna tem- 
perature of the main line of the indicated transition and t,^ is the 
optical depth of the main line. 

It can be seen that the calculated rotational (and thus kinetic) 
temperature depends upon the ratio of antenna temperatures mea- 
sured at the peaks of (1,1) and (2,2) line emission, in addition to 
the optical depth of the (1,1) transition. In order to determine the 
range of for which we may wish to define column density in star- 
forming regions, we have examined the ammonia maps of Perseus 
and W3 presented in ? in addition to other works. 
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Figure 3. The variation of ammonia column density with T^. 

Figure 2. Plane of the equilibrium solutions of 7,- over typically observed 

T* 

ranges of ^t"'"'"' and r,„ 
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Measurement of the optical depth of the ammonia (1,1) tran- 
sition requires detection of the hyperfine satellite quadrupoles (see 
Fig.[D. The faintest of these is the quadrupole at ~23.696 GHz, la- 
belled 'A' in Fig.[T]with an antenna temperature ratio relative to the 
main line of 0.26. As described above, deriving rotation (and hence 
kinetic) temperature from ammonia observations at 1.3 cm relies 
upon detection and measurement of the (2,2) main line. Therefore, 
if we are to determine the optical depth and rotation temperature 
from observations of the (1,1) and (2,2) ammonia lines then we 
require ''.0^ > 0.26. 

' Am(l.l) 

Measurements of the optical depth of ammonia in W3 and 
Perseus indicate that T,n < 5.0 in these regions (?). However, an ex- 
amination of previous observations indicate that it may be as high 
as 10.0 in some star-forming regions (??) (values of may be 
higher than this in actuality but would likely be unobservable due 
to line saturation). The lowest values of previously observed in 
star- forming regions are ~0.1 (??) though it should be noted that 
values of t,^ below ~0.5 become difficult to measure through line 
strength ratios as the ratios approach the theoretical limit. Mod- 
elling and fitting of the hyperfine structure of the ammonia inver- 
sion transitions (e.g. ??) may provide lower values however. 

Solving for over the given ranges of and T,n allows us 

Am(l,l) 

to identify a plane of solutions applicable to star-forming regions. 
The equilibrium solutions of this plane are shown in Fig. |2] and 
represent the range of values of which may be expected to be 
found in any non-extreme star-forming region. T,. is found to vary 
between 8.7 and 32.8 K, corresponding to a kinetic temperature 
range of 8.9 - 45.5 K. 



3 DISCUSSION 

The various ways in which the assumed value of afi^ects the cal- 
culated column density were addressed in section |2] The variation 
of column density over the expected kinetic temperature range de- 
scribed in Section|2](8.9 - 45.5 K) was examined for a grid of values 
of linewidth and optical depth (0.5 km/s < dV < 2.0 km/s, 0.1 < t < 
10.0) in order to investigate the effect of under- or over-estimating 
T^. The resulting change in column density was found to be as high 



as 80% around the midpoint for any given set of dV and t values. 
This uncertainty is relevant to the extreme case; within any given 
region observed values of Tj^ typically vary by a smaller amount, 
e.g. for the Perseus low-mass star-forming region, values of 
range from ~9K to ~20K, while the W3 high-mass star-forming 
region exhibits a range of ~ 12K <Ti^< ~35K (?). The standard de- 
viations of the kinetic temperature distributions in these regions are 
2.3 and 4.9 K for Perseus and W3 respectively. Using these stan- 
dard deviations a more realistic 1 cr variation for the column density 
is -30% 

This uncertainty is equivalent to the measurement error in the 
estimation of optical depth, assuming an antenna temperature mea- 
surement error of 20% which is appropriate for the observations 
used in this work and presented in ?. 

The essentially small variation in the derived column density 
related to temperature means that column density can be calculated 
with the assumption of a mean temperature without any great loss 
of precision. This then allows column density to be determined for 
regions which do not have adequate SNRs in the (2,2) transition for 
kinetic temperature measurements. 

? presented a database compiled from the literature of 264 
star-forming cores observed in the (1,1) and (2,2) lines of ammonia. 
They found a mean of 19 K, in agreement with values found in 
our own data and in other cited works. However, they note that this 
average temperature varies from region to region and attribute this 
variation to association of sources with clusters, where clusters are 
associations of at least 30 embedded stars. Fig. |4] shows the distri- 
butions of values of T^^ derived from our KFPA maps in Perseus and 
W3, along lines of sight where both (1,1) and (2,2) detections were 
obtained. Ranges, averages and standard variations of the T^^ values 
for both regions are listed in Table [T] The plots show an approxi- 
mately Gaussian distribution of values in W3 and an apparently 
skewed distribution in Perseus. It is noteworthy that the mean and 
median averages of the distribution in W3 are significantly higher 
(~17K) than for Perseus (~ 12K). It is possible that the distribution 
of values of Tj^ is truncated in Perseus due to the lower average. 
Below temperatures of ~9 K (2,2) to (1,1) line ratios fall below 0.3 
and thus detection of the (2,2) line becomes increasingly difficult. 
Without a detection of the (2,2) line it is not possible to determine 
and so our sampling of values below 9 K is likely to be highly 
incomplete. 
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Figure 4. The distributions of values of tlie kinetic temperature for maps of 
clumps in Perseus (left) and W3 (right). 



Table 1. Statistics of the Tj^ and N^^j to N^oi Ratio distributions for Perseus 
and W3 clumps. 
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Min 


8.90 


8.90 


0.16 


0.30 


Max 


31.00 


51.50 
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11.01 


Mean 


12.84 


18.73 


1.07 


1.06 


Median 


12.17 


17.66 


1.08 


1.04 


Std.Dev 


2.65 


5.86 


0.25 


0.32 



3.1 Average-Temperature Column Density 

We have estimated ammonia column densities from our KFPA 
data using both the traditional method and the proposed average- 
temperature method. The proposed average-temperature methods 
uses the mean kinetic temperature specific to the detected ammo- 
nia emission in each KFPA map. The results are presented in Fig[5] 
where it can be seen that the area over which column density is de- 
termined increases by 20-30% when using the average-temperature 
method. This allows the column-density maps to extend to lower 
levels, with a minimum column density lower by a factor of ~2 
compared to the traditional determination of column density. The 
expansion of mapped column density, tracing the distribution of 
mass in dense regions, will be beneficial to star-forming studies. 
Accurate determinations of core sizes and tracing the radial mass 
distribution associated with ammonia cores is important in deter- 
mining both virial ratios and the association of gas mass with dust 
mass. By minimizing the lowest known column density contour we 
can improve the quality and accuracy of such studies. 

Where both the (1,1) and (2,2) lines were detected, we can 
compare the results of the two methods directly. Fig |6] shows the 

N' 

distribution of (fj^), where N^^j is the result of the approximate 
method and N(.oi of the exact method, along each line of sight, for 
both Perseus and W3. 

It can be seen that N^^, does not diverge from N^oi by much 
more than 20% over the majority of observed lines of sight and that 
the average-temperature method introduces an additional 1-sigma 
uncertainty of ~ 25 - 30%. Given that the uncertainties in N^oi due 
to measurement error are typically this size (see ?), the approximate 
method does not introduce significant additional error. 
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Figure 6. The distiibutions of values of the average-temperature column 
density (N^.^,]) to traditionally-determined column density (N^oi) ratio. The 
values over mapped clumps in Perseus are shown on the left and values for 
mapped clumps in W3 are shown on the right. 



As the average value of T^^ varies between different star- 
forming regions (Fig. [4] ?) we can minimise the introduced eiTor 
by adopting a mean value and range appropriate to the region 
under study. For example, an average T^^ determined from a cata- 
logue of arbitrarily chosen star-forming cores may be representa- 
tive of physical conditions on a Galactic scale. However, this may 
have little relevance to the column density of a particular region. In 
the presented examples the adopted was the mean over the indi- 
vidual mapped region (typically each KFPA map is circular with a 
radius of 5 arcminutes). In practice, the representative temperature 
which should be used will depend on the data at hand, although it 
has been shown (Section[3} that any introduced uncertainty is less 
than 50% and probably less than 25%. 

Given that the (2,2) transition of ammonia is easily simultane- 
ously observable with modem instrumentation when making (1,1) 
observations, determinations should always be possible for at 
least part of an observed region. It is suggested then that use of the 
proposed method should not be used to replace knowledge of ki- 
netic temperature but that knowledge of the distribution of for 
part of a given region may be used to expand the area over which 
column density is determined and probe a factor of two deeper. 

3.2 Systematic Effects 

It is possible that, by increasing the extent of column density map- 
ping in star-forming regions using the described method, system- 
atic biases on the values of measured column density may be in- 
troduced. There are two main causes of concern, firstly that the 
lower column density values being mapped with the temperature- 
averaging method actually represent lower volume density regions. 
Such regions may have significantly lower excitation temperatures 
due to the volume density being lower than the critical threshold, 
leading to systematic enws. Such an effect would likely be identi- 
fied by a significant decrease in the optical depth of those regions. 
Values of optical depth below t,^(i,i) ~0.5 become difficult to mea- 
sure using the line temperature ratio method described by ? and so 
values less than this can be used as a somewhat arbitrary limit, be- 
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Figure 5. Sections of the Perseus and W3 ammonia maps presented in ? comparing the traditional determination of column density (left or bottom) to the 
described average-temperature determination (right or top). 
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low which associated values of temperature-averaged column den- 
sity might become suspect. An examination of the optical depth 
values of the regions presented here indicate that this limit is not 
violated in any significant portion of the extended maps. 

The other possible concern is that maps of column density us- 
ing the proposed temperature- averaged method might be extended 
to the point that cores may be externally heated, by PDR heating 
in the external portions of the clouds for example. This effect has 
been observed by ? in measuring the excitation temperatures of 
CO isotopologues at varying optical depths in star-forming regions 
close to HII regions. If this were to occur, then the assumption of 
an average temperature associated with a star-forming core would 
become invalid, again leading to systematic errors in the derivation 
of column density. 

External heating is unlikely to be an issue in the determination 
of temperature-averaged column density for two reasons; firstly, it 
can be seen from Fig. |5] that the overall expansion of regions re- 
sulting from the assumption of an average temperature is small in 
comparision to the overall size of the regions. Thus, any effects of 
external heating would be evident in the observed values of T-^ and 
column-density determined in the usual manner. Secondly, the re- 
sults of ? show that C"*0 temperatures are relatively independent 
of external heating and thus internal, super-critical density regions 
associated with low optical depth tracers are likely shielded or oth- 
erwise unafi"ected by significant external heating. 



4 CONCLUSION 

The proposed method of deriving column density through the 
assumption of an average or representative temperature provides a 
way for the mass distribution in star formation regions to be traced 
over a significantly greater regime than previously possible. While 
this method cannot be expected to be a substitute for knowledge of 
the kinetic temperature, the ever-growing prevalence of, and need 
for, fast mapping of large regions means that it can be a useful tool. 
The accuracy of column density measurements made through the 
assumption of a mean temperature has been shown to be compa- 
rable to measurements made with full knowledge of temperature. 
Also, the additional uncertainty introduced is likely to be less than 
that from measurement errors. These extended regions of emission 
are important to studies of star formation in many ways, not least 
in studies of the relation of submillimetre continuum emission to 
gas mass and the distribution of virial ratios in different regions (?). 
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